Multiband superconductivity in PrPt 4 Gei 2 single crystals 
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We report measurements of the London penetration depth AA(T) and the electronic specific heat 
C e (T) on high-quality single crystals of the filled-skutterudite superconductor PrPt4Gei2 (T c ~8K). 
Both quantities show a weak temperature dependence at T <C T c , following AA ~ T" (n ~ 3.2) 
and C e /T ~ T 2 ' 8 . Such temperature dependences deviate from both conventional s-wave type and 
nodal superconductivity. A detailed analysis indicates that the superfiuid density p s {T), derived 
from the penetration depth, as well as the electronic specific heat can be consistently described in 
terms of a two-gap model, providing strong evidence of multiband superconductivity for PrPt4Gei2. 

PACS numbers: 74.25.Bt; 74.20.Rp; 74.70.Dd; 74.70.Tx 
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The filled-skutterudite compounds MT^X^ (M=rare- 
earth or alkaline-earth metals, T=Fe, Ru, Os, and X~P, 
As, Sb) demonstrate remarkably rich physical proper- 
ties [l(. Particular attention has been paid to super- 
conductivity (SC) observed in the Pr-based compounds. 
For example, PrOs 4 Sbi2 is a heavy-fermion superconduc- 
tor with T c — 1.85 K [2j. Electrical-quadrupole, rather 
than magnetic-dipole, fluctuations are believed to medi- 
ate the Cooper pairs in this compound, which is unique 
among heavy fermion superconductors. The supercon- 
ducting order parameter of PrOs 4 Sbi2 remains contro- 
versial: nodal SC 0, 0] as well as s-wave SC [&] were 
proposed. More recent experiments seem to support a 
scenario of multiband SC [6j, lZ|- On the other hand, the 
isostructural compounds PrRu4Sbi2 and PrRu4Asi2 ap- 
pear to be s-wave superconductors [1, Q . 

Recently, a series of new skutterudite superconductors 
with a germanium-platinum framework, i.e., MPt4Gei2 
(M=Sr, Ba, La, Pr), were successfully synthesized 
[HA Among all the Pr-filled variants, PrPt4Gei2 

shows an unexpectedly high transition temperature of 
T c =7.9 K [IH. The Sommerfeld coefficient of PrPt 4 Gei 2 
(7„ = 76 mJ/mol K 2 ) [l2| is comparable to that of 
PrRu4Sbi2 [U and PrRu4Asi2 @, but much smaller than 
that of PrOs4Sbi2 [2i- Furthermore, the crystalline elec- 
tric field (CEF) splitting of the J = 2 Hund's rule mul- 
tiplet of Pr 3+ between the ground state and the first 
excited state is rather different among these Pr-based 
superconductors, e.g., Aqef = 7 K for PrOs4Sbi2 [HI 
and Acef =130 K in PrPt 4 Gei 2 (HQ. It is, there- 
fore, of great interest to systematically compare the su- 
perconducting properties of these materials, which may 
help to elucidate their pairing mechanisms. Similar to 
PrOs4Sbi2, previous studies on polycrystalline samples of 
PrPt 4 Gei2 showed controversial results. Measurements 
of the specific heat and muon-spin rotation (/iSR) sug- 
gest the possible existence of point nodes in the super- 
conducting gap [12] ; zero- field ^SR also provides evidence 



of time reversal symmetry breaking below T c |l5| . simi- 
lar to what was observed for PrOs4Sbi2 jlH- However, 
73 Ge nuclear quadrupole resonance (NQR) experiments 
display a pronounced coherence peak in the spin-lattice 
relaxation rate 1/Ti at temperatures just below T c , sug- 
gesting s-wave SC [l7| . Very recently, a possible scenario 
of multiband SC was proposed for PrPt4Gei2, based on 
the analysis of the critical fields [l8[ as well as photoe- 
mission spectroscopy [ijjj]. However, these experiments 
were performed on polycrystalline samples at relatively 
high temperatures, which could not make a clear asser- 
tion on the gap symmetry. The reasons underlying such 
discrepancies of the gap structure in PrPt4Gei2 are not 
yet clear, and further measurements, in particular those 
based on high-quality single crystals, are badly needed. 

In this Letter, we probe the superconducting gap sym- 
metry of PrPt 4 Gei2 by measuring the London penetra- 
tion depth AA(T) and the specific heat C P (T) of high- 
quality single crystals. Precise measurements of the 
penetration depth changes at low temperatures show 
AA ~ T n with n ~ 3.2, indicating that PrPt4Gei2 is ac- 
tually neither a simple BCS nor a nodal superconductor. 
A detailed analysis of the superfiuid density p s (T), con- 
verted from A(T), and the electronic specific heat C e (T) 
provide strong evidence of two-band SC for PrPt4Gei2- 

High-quality single crystals of PrPt4Gei2 were syn - 
thesized by using multi-step thermal treatments [20j. 
Powder X-ray diffraction indicates the presence of a 
small amount of foreign phases. Energy-dispersive X- 
ray (EDX) analysis confirms that all the crystals have 
a stoichiometric composition and the impurity phases, 
mainly PtGe2 and free Ge, are located at the crystal 
surfaces (20| . In our measurements, the crystals were 
mechanically polished to get rid of these surface con- 
taminations. Precise measurements of the resonant fre- 
quency shift A/(T) were performed by utilizing a tunnel 
diode oscillator (TDO) based, self-inductance method at 
an operating frequency of 7 MHz down to about 0.5K 
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in a 3 He cryostat [2l|. The change of the penetration 
depth is proportional to the resonant frequency shift, 
i.e., AA(T) = GA/(T), where G is solely determined 
by the sample and coil geometries In this context, 
AA(T) is extrapolated to zero at T = by polynomial 
regression, i.e., AA(T) = A(T) — A . Here the value 
of zero-temperature penetration depth, Ao=114 nm, was 
adopted from previous /iSR experiments [12j . Measure- 
ments of the magnetic susceptibility were carried out in a 
SQUID magnetometer (Quantum Design), and the heat 
capacity was measured in a 3 He cryostat, using a relax- 
ation method. 

Figure 1 presents the electrical resistivity p(T) for 
PrPt4Gei2, which shows an S-shape behavior upon cool- 
ing down from room temperature, as often observed in 
d-band materials. A sharp superconducting transition, 
evidenced from both the electrical resistivity p(T) (inset 
(a)) and the magnetic susceptibility x(T) (inset (b)), to- 
gether with a large resistivity ratio (/9(300K)/p(8K)=19) 
confirm the high quality of our single crystals. Further- 
more, the superconducting transition temperatures T c , 
determined from the zero resistivity and the onset of the 
magnetic susceptibility, are nearly the same (T c ~7.8 K), 
proving good homogeneity of the samples. 

The inset of Fig. 2 shows the change of the penetration 
depth AA(T) up to 9K, which reveals a sharp supercon- 
ducting transition at T c ~ 8.1 K. Here G=2.13A/Hz. It 
is noted that we have measured the penetration depth 
for several samples and the data are well reproducible. 
The values of T c , derived for different samples by dis- 
tinct methods, are nearly same too. In the main plot of 
Fig. 2, we presents AA(T) at low temperatures, together 
with the fits of various models to the data. Obviously, 
the standard BCS model can not describe the experi- 



mental data. Moreover, the penetration depth AA(T) 
deviates also from that of nodal SC, for which a linear 
and quadratic temperature dependence is expected for 
the case of line and point nodes, respectively. Instead, a 
power law of AA ~ T n (n ~ 3.2) presents a reasonable 
fit to the experimental data. An enhanced power-law ex- 
ponent n, e.g., a quadratic temperature dependence in 
d-wave superconductors, may arise from nonlocal effects 
or impurity scattering [22j. However, such possibilities 



are excluded for PrPt4Gei2 because both the penetra- 
tion depth (Ao =114 nm) 12j and the mean free path 
(I = 103 nm) are much larger than the coherence length 
(£o = 13-5 nm) 12|, implying that the samples are in 
the clean and local limit. Here we estimate the mean 
free path from I = [ \ 8 ~ ^ * ^ I - ) 1 / 2 !23|, where p , 
£o and 7„ represent the electrical resistivity at T C =7.8K 
(po = 3.5 x 1CP 6 Hem), the aforementioned coherence 
length and the Sommerfeld coefficient at T c (j n =1795 
erg cm _3 K -2 ), respectively. On the other hand, multi- 
band effects may also give rise to power-law-like behavior 
at low temperatures with a large exponent n, which will 
be further elucidated by the analysis of both the super- 
fluid density and specific heat. 

The superfluid density p s {T) provides an important 
characterization of the superconducting gap symmetry. 
Fig. 3 shows the temperature dependence of the nor- 
malized superfluid density p s {T) for PrPt4Gei2 (circles), 
which is calculated by p s = [Ao/A(T)] 2 . For compari- 
son, in Fig. 3 we also include the superfluid density from 
the /iSR results determined on polycrystalline samples 
(diamonds) [12]. Obviously, these two data sets are 
quite compatible, although the /iSR data are more scat- 
tered, with a poor resolution when compared with the 
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FIG. 1: Temperature dependence of the electrical resistivity 
p(T) for PrPt4Gei2. Insets show the superconducting transi- 
tions in the electrical resistivity p(T) (a) and magnetic sus- 
ceptibility x(T) (b), respectively. 
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FIG. 2: Temperature dependence of the penetration depth 
AA(T) for PrPt4Gei2. Solid, dash-dotted and dashed lines 
represent the fit of AA ~ T 3 ' 2 , AA ~ T 2 (point node) and 
single-gap BCS model, respectively. Inset shows AA(T) over 
the entire temperature range. 
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TDO results. This allows to probe the gap structure of 
PrPt4Gei2 in a much more precise way than before. 
The superfluid density can be calculated by: 



p s (T) = 1 + 2- 



df 



E 



dE _ A 2 (T) 



dE) 



FS • 



(1) 



where / = (eV £2 + A l( T )/ fc « T + l)- 1 i s the Fermi distri- 
bution function and (. . .)fs denotes the average over the 
Fermi surface. For the temperature dependence of the 
energy gap, we take A(T) = A tanh ^ bT " [2 ACc < T 

I)] ' 5 



24] 



. r2 AC. i T 

A L3 ~f n T a \T c 

is the specific heat jump at T c . 



Here AC e 

Note that Eq. 1 is applicable for various gap functions 
Afe (=A(6, 4>)) in the pure/local limit. Given a gap func- 
tion A(0, </>), then one can fit it to the experimental data. 
Here 9 and <f> denote the angles away from the z-axis 
and x-axis in k-space, respectively. In this analysis, the 
zero-temperature gap amplitude, Ao, is the sole fitting 
parameter. 

Possible symmetries of the order parameter have been 
theoretically investigated for the skutterudite supercon- 
ductors with tetrahedral point group symmetry (Th) pB| . 
Various gap functions A(0, 0), restrained by the crystal 
symmetry, have been adopted to fit the superfluid den- 
sity p s {T) of PrOs4Sb 12 [4]. In this context, we apply 
a similar analysis to the experimentally obtained p s (T) 
data of PrPt4Gei2- Fig. 3 presents the fits of different 
gap functions allowed by the crystal symmetry; the de- 
rived fitting parameters of Ao are summarized in Table 
1. Apparently, the gap functions B and E cannot re- 
produce the experimental data. On the other hand, the 
fits of functions A, C and D are close to the experimental 
data, but show significant deviations at low temperatures 
(inset of Fig. 3). We note that models C and D, both hav- 




FIG. 3: Superfluid density p s (T) versus normalized tempera- 
ture T/T c . Inset expands the low temperature region. Circle 
(o) and diamond (©) dispaly the experimental data derived 
from TDO measurements (this study) and ^tSR experiments 
(from Ref. [12]), respectively. Lines shows theoretical fits of 
various gap functions as listed in Table 1. 



TABLE I: Summary of various gap functions and so derived 
fitting parameters Aq. 



model 


Gap function A(9,<p) 
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1.76 
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| Ao sin 6 sin <j>\ 
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A sin0| 


2.7 
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A (l - sin 4 (9 cos 4 (9) 


2.9 
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A (l - (sin 4 4> + cos 4 <j>) sin 4 6 


) 3.2 



ing point nodes in the superconducting gap, were previ- 
ously assumed to present a good fit to the pSR data [ll| , 
which are rather scattered at low temperature. The more 
precise measurements of the penetration depth, AA(T), 
and the corresponding superfluid density, p s {T), indicate 
that the conventional one-gap BCS model as well as the 
nodal-gap model D provide a poor fit to the low tempera- 
ture data. The nodal-gap model C can fit the TDO data 
relatively well, but significant deviations remain below 
0.4T/T C . Instead, the two-gap BCS model gives the best 
fit to our experimental data. 

In the case of two-gap BCS superconductors, the su- 
perfluid density can be extended to the following linear 
combination[26j : 

p s {T) = xp s (Al,T) + (1 - x)p s (Aq, T), (2) 

where A l (i — 1,2) represent the size of two isotropic 
gaps at zero temperature, and x is the relative weight 
of the contributions from Aj. As shown in Fig. 3, the 
two-gap BCS model nicely fits the experimental data of 
PrPt4Gei2 over the entire temperature region. The so 



derived parameters of Ai = 0.8fcsT c 



A 2 



2.0fc B T c and 



a;=0.15 meet the theoretical constraints that one gap is 
larger than the BCS value and the other one is smaller 
27j . as demonstrated in the prototype two-gap BCS su- 
perconductor MgB 2 



26|. 



In Fig. 4, we present the low-temperature specific heat, 
C P (T), of a PrPt4Gei2 single crystal. A sharp supercon- 
ducting transition is observed at T c = 7.7 K, being close 
to that determined from other experiments. The spe- 
cific heat data above T c can be fitted by a polynomial 
expansion C P (T) = 7 „T + /3T 3 . Here C e = 7„T and 
C p h — /3T 3 denote the electronic and phonon contribu- 
tions, respectively. This yields the Sommerfeld coefficient 
in the normal state, j n = 69 mJ/mol K 2 , and the Debye 
temperature Qd =190 K, which are close to those found 
in case of polycrystalline samples [l2j]. For polycrystals 
12| , a pronounced upturn was previously reported in the 
low-temperature specific heat C e (T)/T. Similar specific 
heat anomalies were also observed in some as-grown sin- 
gle crystals. A careful examination showed that such 
an upturn in C e (T)/T has to be attributed to a nuclear 
Schottky anomaly caused by the Pr-containing surface 
contaminations 20]. Indeed, the specific heat anomaly 
disappears for the polished single crystal as shown in 
Fig. 4(b), allowing us to accurately analyze its low tern- 
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perature behavior. 

The electronic specific heat of PrPt4Gei2, obtained af- 
ter subtracting the phonon contributions, is presented in 
the main part and inset (b) of Fig. 4 as C e /T verse 
T, together with the fits of various models. As shown 
in the inset (b), the data can be well described by a 
power law, C e /T ~ T 28 . This behavior deviates from 
the quadratic temperature dependence of C e (T)/T re- 
ported in Ref. 12]. The discrepancy is likely to result 
from the nuclear Schottky anomaly of the polycrystalline 
samples discussed before. With the previous data [Hj], a 
proper subtraction of this Schottky anomaly is difficult 
and, therefore, deviations from the true specific-heat be- 
havior become likely at low temperatures. Furthermore, 
the standard BCS-model is not sufficient to fit the exper- 
imental data (Fig. 4(b)), while the two-gap BCS model 
presents the best fit to the C e (T)/T data (main figure). 
According to the phenomenological two-gap BCS model, 
the heat capacity is taken as the sum of contributions 
from the two bands, each one following the BCS-type 
temperature dependence [28| . In the main panel of Fig. 
4, we plot the contributions from the two superconduct- 
ing gaps, Aj — 0.8fcsT c and A 2 , = 2.0fcsT' c , as well as 
their sum (solid line). The weight contributed from the 
first gap, Aj, is x=0.12. All these fitting parameters are 
remarkably consistent with those obtained from the su- 
pcrfluid density p s (T) , providing strong evidence of two- 
gap SC for PrPt 4 Gei 2 . 

Evidence of BCS-like SC, including two-gap type, has 
been observed in several skutterudite compounds. For 
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FIG. 4: Temperature dependence of the electronic specific 
heat C e (T)/T of a PrPt4Gei2 single crystal at zero field. 
Dash dotted lines (green and blue) and solid line show the 
individual and total contributions of the two gaps to the 
specific heat C e (T)/T, respectively. Inset (a): total specific 
heat, C P (T)/T, plotted as a function of T 2 . Inset (b): elec- 
tronic specific heat C e (T)/T at low temperatures, fitted by 
C e /T ~ T 2 ' 8 (dotted line). Dashed line refers to the standard 
BCS model. 



example, PrRu4Sbi2 0, PrRu4Asi2 Q and their non- 
47 counterparts [f2[ are believed to be s-wave supercon- 
ductors. Recent measurements indicate that PrOs4Sbi2 
is an extreme two-band superconductor [f| 0] ; here en- 
ergy nodes were assumed to exist in the small gap, and 
the isotropic large gap dominates the superconducting 
properties near T c [7|, or when a sufficiently large mag- 
netic field is applied (29| . Tw<>gap BCS SC was also 
proposed for both PrRu4Sbi2 0| and LaOs4Sbi2 [30j . 
the latter one suggesting that 4/-electrons are not the 
origin of multiband SC. In PrPt4Gei2, band structure 
calculations indicate an only minor contribution of the 
4 /-electrons to the density of states at the Fermi energy, 
suggesting that the 4/-electrons may not be playing a 
significant role on SC in this compound either In- 
deed, the thermodynamic properties and low-lying CEF 
scheme of PrPt4Gei2 seem to be rather different from 
those of the heavy fermion compound PrOs 4 Sbi2, but re- 
semble other skutterudite compounds [H 0] . Indications 
of two-gap SC for PrPt4Gei2 were also inferred from re- 
cent measurements of the upper and lower critical fields 
(ljj and photoemission spectroscopy [l9|. Furthermore, 
multiband SC is compatible with the observations of a 
coherence peak in the NQR measurements just below 
T c 17j. Such a multi-gap structure seems to be char- 



acteristic for the skutterudite superconductors; the small 
gap, either with or without nodes, is rather subtle and 
can be easily destroyed by external effects, e.g., a mag- 
netic field, so that the large gap is predominant. Recent 
/iSR measurements performed on polycrystalline samples 
of PrPt4Gei2 showed evidence of time-reversal symme- 
try breaking 15j. To confirm it and check the possi- 
ble existence of nodes in the small gap of PrPt4Gei2, it 
would be important to repeat these measurements with 
high-quality single crystals. Detailed calculations of its 
electronic structure are also highly desirable in order to 
further elucidate the multiband structure in PrPt4Gei2. 
Moreover, comparative studies of the Pr-based skutteru- 
dites and the non-/ electron isostructural compounds, 
e.g., MPt4Gei2 (M = Sr, Ba and La), are necessary to 
reveal the potential role of /-electrons on SC. 

In summary, we have studied the superconducting or- 
der parameter of PrPt4Gei2 by measuring the penetra- 
tion depth AA(T) and specific heat C P (T) on high-quality 
single crystals. For T <C T c , both quantities demon- 
strate a weak temperature dependence and can be fit- 
ted by a power- law behavior with a large exponent, i.e., 
AA ~ T 3 - 2 and C e /T ~ T 2 8 , which is inconsistent with 
both a single-gap BCS model and nodal-gap SC. Instead, 
we can describe the superfluid density p s (T) and the elec- 
tronic specific heat C e (T) in terms of a phenomenologi- 
cal two-gap BCS model with consistent gap parameters 
of AJ = 0.8k B T c , A 2 = 2.0k B T c and x = 0.12 - 0.15, 
the weight contributed by the small gap. These findings 
have elucidated the controversial results found in the lit- 
erature and provide unambiguous evidence of multiband 
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SC for PrPt 4 Gei 2 . 
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